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Abstract: We investigate the lepton flavor violation (LFV) in the inert scalar 
model with higher representations. We generalize the inert doublet model with right 
handed neutrino by using higher scalar and fermion representation of SU{2)i. As the 
generalized model and the inert doublet model have the same parameter space, we 
compare the rates of /r —)■ ey, fx —)■ eee and fx — e conversion in nuclei in the doublet 
and its immediate extension, the quartet model. We show that the corresponding 
rates are larger in the case of higher representation compared to the Inert doublet 
for the same region of parameter space. This implies that such extended models are 
more constrained by current LFV bounds and will have better prospects in future 
experiments. 
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1 Introduction 

Neutrino oscillation provides the direct evidence for lepton flavor violation in the 
neutrino sector. Therefore, one also expects LFV in the charged lepton sector which 
is yet to be observed. This is a generic prediction in most of the neutrino mass mod¬ 
els and depending on the realization details of the model, the rates of different LFV 
processes can be very different. In this paper, we have focused on radiative neutrino 
mass model at one loop proposed in [1], known as the scotogenic model, where the 
scalar content of the model is the inert doublet. Apart from its role in neutrino mass 
generation, the inert doublet has been extensively studied in the context of dark 
matter [2-9], mirror model and extra generation [10, 11], electroweak phase transi¬ 
tion [12-16] and collider studies [17, 18, 21]. As the higher scalar representation is 
not forbidden by any symmetry in the model, the immediate generalization of the 
doublet, the quartet with isospin J = 3/2 was studied in [22] to check whether it is 
viable in providing both light scalar dark matter and strong electroweak phase tran¬ 
sition in the universe. Here we have incorporated higher scalar representation instead 
of the doublet in the scotogenic model and determined the viable SU{2)l fermion 
multiplet for generating neutrino mass. LFV processes in the scotogenic model with 
inert doublet has been studied in [23-28] (and references therein). The extension 
of the scotogenic model has been addressed in [29, 30]. Also larger multiplets have 
been incorporated in type III seesaw model [31] and in models of radiative neutrino 
mass generation at higher order with dark matter [32]. 

The generalization of scotogenic model with higher SU{2)i half-integer repre¬ 
sentation does not change the parameter set of the Lagrangian of the inert doublet 
at the renormalizable level. Therefore it gives us the opportunity to investigate the 
predictions of LFV processes for different scalar representations for the same region 
of parameter space. In particular, we compare the LFV processes for the doublet 
and the quartet in the light of current experimental bounds and future sensitivities. 

There have been many great experimental efforts to detect positive LFV signal in 
la —t IjS'j, la —t 3//3 and fi — e conversion rate in nuclei. In the case of muon radiative 
decay, the MEG collaboration [33] has put a limit of Br(/r —>■ ey) < 5.7 x 10“^^ 

[34] and will have sensitivity of 6 x 10“^^ after acquiring data for three more years 

[35] . In addition, current bound on branching ratio of lepton flavor violating 3-body 
decay, /i —)■ eee is 1 x 10“^^ set by SINDRUM experiment [36] and Mu3e experiment 
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will reach a sensitivity of 10“^® [37]. Furthermore, SINDRUM II experiment has put 
current limit on muon to electron (/i — e) conversion rate in Gold (An) and Titanium 
(Ti) nucleus of 7 x 10“^^ [38] and 4.3 x 10“^^ [39] respectively. The future projects 
Mu2e [40, 41], DeeMe [42], COMET [43] and PRISM/PRIME [44, 45] will improve 
this bound from 10“^^ to 10“^®. For other LFV processes and their experimental 
bounds, please see Table I of [27]. We have compared the predictions of the LFV 
processes fi —)■ ey, fi —>■ eee and p — e conversion rate in An and Ti for both doublet 
and quartet scalars and our comparison has revealed that the contributions of the 
quartet in all LFV processes are larger than those of the doublet for the same region 
of parameter space. Consequently, the contribution of higher scalar representation 
to LFV processes have better experimental prospects. 

The paper is organized as follows. We describe the model in section 2. In section 
3 we present the relevant formulas of /i —)■ ey, fi —)■ eee and fi — e conversion processes 
for the inert doubler and quartet. We present the result in section 4 and conclude in 
section 5. appendix A contains the mass spectrum of the inert doublet and quartet 
in our parametrization. The expressions of the loop functions are given in appendix 
B. In appendix C we collect the Feynman diagrams for /rey vertices, fieZ vertices 
and box diagrams. 

2 The Model 

Any multiplet charged under SU{2)l x U{1)y gauge group is characterized by the 
quantum numbers J and V, with the electric charge of a component in the multiplet 
is given by Q = T 3 + V. For half-integer representation J = n/2, T 3 ranges from 
— I to |. So the hypercharge of the multiplet needs to be V = ±T 3 for one of the 
components to have neutral charge. For integer representation n, similar condition 
holds for hypercharge. 

The generalized scotogenic model involves one half-integer SU{2)l scalar multi¬ 
plet A with hypercharge V = 1/2 and three generations of real {Y = 0) odd dimen¬ 
sional fermionic multiplets, Fj (i = 1 — 3) charged under Z 2 symmetry, A —)■ —A 
and Fi —)■ —Fj. When the scalar multiplet is hxed to be J = n/2 , n odd, there are 
two choices for fermionic multiplet which can give Z2 even SU{2)l x U{1)y invariant 
Yukawa term with the lepton doublet; J = or The charged lepton sector is 
augmented by the following terms 

Mp. _ — 

C D- ^F^P^Fi + yiaFi-la-^ + h.c (2.1) 

where the dot represents the proper contractions among SU{2) indices. In the sub¬ 
sequent analysis we have chosen fermion multiplet to be J = 
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The general Higgs-scalar multiplet potential , symmetric under Z 2 , can be writ¬ 
ten in the following form, 

Ho($, A) = + MqWA + + A2(A'fA)2 + AslA^T^Ap + a$t$AtA 

+ /3$V“<hA^T“A + 7 [($'^er“$)(A^C'T“A)'^ + h.c] ( 2 . 2 ) 


Here, r“ and T“ are the SU{2) generators in fundamental and A’s representation 
respectively. C is an antisymmetric matrix analogous to charge conjugation matrix 
dehned as, 

CT^C-^ = -T^t ^ 2 . 3 ) 

Since C, is an antisymmetric matrix, it can only be dehned for even dimensional 
space, i.e only for half-integer representation. If the isospin of the representation is J 
then C is (2 J -|- 1) x (2 J -|- 1) dimensional matrix. The generators are normalized in 
such a way so that they satisfy, for fundamental representation, Tr[r“r^] = and 
for other representations, Tr(T“T^) = iA 2 (A) 5 “^. Also T“T“ = C 2 (A). Here, T* 2 (A) 
and 6*2 (A) are Dynkin index and second Casimir invariant for A’s representation. 
Notice that, 7 term is only allowed for representation with {J,Y) = (|, |) and it is 
essential for the generation of neutrino mass at one-loop. 

The scalar representation with {J,Y) = (|, |) and the fermionic representation 
with {J,Y) = (^^,0) have the component helds denoted as A*^*5) ^^d F^Y) respec¬ 
tively where Q is the electric charge. They are written explicitly as 



/ A(^) \ 


/fy-^^\ 

An = 

2 

A(o) = j^(S + iA) 

and F n-i = 

2 

_p(o) 


A(-V) J 




For the former representation every component represents a unique held while for 
the latter there is a redundancy fYQ) = (^F^'^^)*. 

The choices for real fermion multiplet with the doublet are either (J, Y) = (0, 0) 
or (1, 0) and with the quartet, choices are either (J, Y) = (1, 0) or (2,0). Our analysis 
has focused on the following pairs of scalar and fermionic multiplets: (Aj^i,Fij=o) 
and (Aj^ 3 ,In component helds, the doublet scalar D, right handed (RH) 
neutrino, and the quartet scalar A and the triplet fermion Fj are expressed as 


D 






I \ 

A+ 

AO ^^(^ + ,71) 

V / 



i 


(2.5) 
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2.1 Mass spectra 


We now sketch the general form of mass spectrum for the scalar and fermionic multi- 
plet which was also presented in [22] . The neutral component of the scalar multiplet 
{Y = 1/2) will have T 3 eigenvalue as T 3 = — 1. Now for the Higgs vacuum expecta¬ 
tion value, (<h) = (0, the term (<hl)r^(<h) gives —So masses for the neutral 
components, S and A are splitted by the 7 term as 

ml = + ^(^a+^[3 + ( 2 - 6 ) 

ml = + ^(^a+^/3 - ( 2 - 7 ) 

Here, p = lDim(|) = 1, 2,... comes from 2p x2p C matrix. For the charged compo¬ 
nent, with T 3 = m, where, m = n/ 2 ,n /2 — 1,..., —n/2, we have, 

rnfm) = + \ (^oi - m^ v\ ( 2 . 8 ) 


Moreover, because of the 7 term, there will be mixing between components 
carrying same amount of charge. A component of the multiplet is denoted as |J, T3). 
Components with ||,m) and ||,—(m -|- 1)) (such that —m — 1 > —|) will have 
positive and negative charge Q = m + ^ respectively. Now (<h)^er“(<F) gives 
Therefore, the mixing matrix between components with charge IQI is. 


( 




m 


(m) 


Vw\/(f-"') (f+ "^ + l) 

And the mass eigenstates are. 


■yv^ 

~T~ 



t-m)(f + m+l)^ 


A'« = cosOq a 7, + siti»Q a7„_i, 

A? = -sm«QA7, + cos«eA7„,_i, 


(2.9) 


( 2 . 10 ) 


where we have 


tan26*Q = 


2 (MJ)i 2 


(A7)„ - (a7)22 


( 2 . 11 ) 


Note that the real fermionic multiplet is degenerate at the tree level. However, 
there is a small splitting between the charged and neutral component due to radiative 
correction which is O(lOOMeV) [46]. This splitting is needed in order to treat the 
neutral fermion as the dark matter candidate. 
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Figure 1. Neutrino mass generation in the inert doublet (first figure from the left) and 
the quartet (second and third figures). 


2.2 Neutrino mass generation 


The light neutrino masses are generated at one-loop level as shown in hgure 1. The 
neutrino mass matrix is expressed as 


m. 


u ja/S 


E 

2=1 


yaiyipMp. 


<-'1 0 -1 
2 2 


mt 


2 2 
Trig — mp_ 


In^ 


m: 


m. 


m\ - m% 


In^ 


rrii 




Qt^o 


m 


Q,i 


Tu¬ 


rn 


Q,i 


m 


Q,2 


-In- 


m: 


Q,2 


m: 


Q,i 


m. 


ml 


m, 


Q ,2 


mt 


m. 


= {y'^^y)afi 


( 2 . 12 ) 


Here Cmi,m 2 ,m 3 is the Clebsh-Gordon (CG) coefficient and mi, m 2 and m 3 are the 
T 3 eigenvalues of lepton doublet, fermion and scalar multiplet respectively. Non-zero 
CG coefficient requires mi -|- m 3 = m 2 . Also Ri^m is the element of the rotation 
matrix that mixes the two scalar components with same charge |( 2 | and m^ ^ are the 
corresponding mass eigenvalues. Moreover, Aj is the loop function. 


A,- = 


Mp. 

IGtt^ 


G?o_r 

2 2 


mt 


-In 


mt 


m. 


-In 


m\ 


mi. 


ml 


ml 


ml 


ml 


ml 


E 




Rl,mR2,—m—l 


m 


Q ,1 


-In- 


m 


Q,i 


m: 


Q ,2 


-In- 


m 


<9,2 


m 


<9,1 


ml 


ml 


m 


< 9,2 


ml 


ml 


(2.13) 


Therefore the neutrino mass at one loop in the doublet case is given by 


im. 


\ doublet 

lap 


E 

2=1 


yaiyif}-^Ni 


ml 


m| - m% 


In 


ml 


m\ 


m 


Ni 


ml, — m 


-In 


ml 


Ni 


m 


Ni 


(2.14) 


where Mtv^ is the mass of the i-th right handed neutrino. When m| ~ m\ = m^ 
then eq. (2.14) gets simplihed 


m. 


\ doublet 

lap 


E 

2=1 


yaiyiyiv 

IGtt^Mtv 




+ 




Tu¬ 


rn 


N,, 


ml 


m 


Ni 


ml 


m 


Ni 


ml 


(2.16) 
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On the other hand, the nentrino mass matrix in the qnartet case is given by 


( A/T 


2 

'y ^ yai^iVip 

i=l 


(2.16) 


with the loop factor, 

y^quartet _ ^ 


3(47r)2 

1 

6(47r)2 


M, 


Fi 


ml. 


mt 




In 


sin 29Mpi 


m 




m 


m. 


K. 


m, 


K. 


In 




m 




Mk 


In 


4- 


m 




Ki 


m 


m, 




Mk 


In 


M'k 


(2.17) 


Explicit expressions of masses in the inert doublet and quartet models are included 
in appendix A. 

The neutrino mass matrix can be diagonalized as 


UpMNS^J'UpMNS = '^v (2-18) 


where 


1 

( C12C13 

5i2Ci3 

5136*^^ 

; j 

^10 0 \ 

UpMNS = 

-S12C23 - Ci2S23Si3e“*'^ 

C12C23 - 5125235136"*'^ 

523C13 

X 

0 0 

\ 

V 512523 - Ci2C23Si3e“*'^ 

-612523 - 5126235136"*'^ 

C23C13/ 

1 1 

{0 0 eO/2/ 


(2.19) 


Here, Qj = cosOij, Sij = sin^jj, 6 is the Dirac phase and a, (3 are the Majorana 
phases. 

The Yukawa matrix Ui^ {a = e, p, r) is expressed using the Casas-Ibarra parametriza- 
tion [49] so that the chosen parameter space automatically satishes the low energy 
neutrino parameters, 

y = \/A R ^/fhy Upj^j^g ( 2 . 20 ) 

where i? is a complex orthogonal matrix. 


2.3 Perturbativity 

If there are N generations of right handed fermion multiplet, perturbativity of the 
Yukawa gives the following constraint [50, 51] 

3 Af 

= Y Y ~ 0(1) (2.21) 

i=l j=l ^ 

If R is taken to be real, the constraint translates into the largest ratio, ^ 0(1)) 
whereas for the general case when R is complex, each entry will be bounded as 
I R. .1 < / Y 

puj I ~ y 3Nrh^. ■ 
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3 Lepton flavor violating processes 


In this section we have presented the relevant analytical formulas of LFV processes 
for the doublet and quartet case. In the standard model due to the GIM suppression 
the rate of /i —)■ ey becomes ~ 10“®^ thus negligible. On the other hand the presence 
of heavy right handed neutrino that mixes with left handed (LH) neutrinos, spoils the 
GIM suppression and one could obtain the rate which can be probed by experiment 
[52-58] . In inert scalar models, Z 2 symmetry forbids the mixing between LH and RH 
neutrinos but the enhancements in the LFV processes are provided by the — N^. 
loops in the doublet and A — Fi loops in the quartet model. We have focused on 
three LFV processes: /i —)■ ey, fi —)■ eee and fi — e conversion in nuclei in this paper 
as they have the most stringent limits from the experiments. 


( 3 , 1 ) 


3.1 p —)■ ey 

The branching ratio for fi —)■ ey, normalized by Br(/i —)■ elyz/^), is [27, 59] 

Br(/r ey) = Br(/i ez/^Iy) 

where Hd is the dipole form factor. The Feynman diagrams of one-loop contributions 
by the doublet and quartet to the /rey vertex that enters into the dipole form factor 
calculation, are given in hgure 9. 

The contributions from the doublet is the following, 

3 * 1 

/(doublet _ \ ^ VeiVitJ- ^ Tr(n) f ^ \ /o n\ 


327r2 my 

i=l ^ 


Here [x) is the loop function given in the appendix B and Xi^ = where 

a = C^. On the other hand, the quartet contribution will have two parts 

^quartet ^quartet , ^quartet /o q\ 

— ^D{n) ^D{c) 

where is the contribution of the neutral component and is that of the 

charged component of the fermion triplet. Also, for the notational convenience, we 
introduce generalized Yukawa coupling yi^a = ViaCa where is the corresponding 
Glebsh Gordon coefficient associated with a-th component of the quartet. The two 
contributions are 

VeiaVif^cr 1 


^quartet _ \ ^ \ ^ 

^Din) - 2^ 

i=l (7 


327r2 m'i 




( 3 . 4 ) 


where Xi^ = a = A^, A^. And 


/■quartet \ ^ \ ^ VeicrVilJ-fr 1 t7'(c)/' 

^Dic) y 

i=l a ^ 

where Xi^ = and a = A++, S, A. 


Xi, 


(3.5) 






3.2 fi —)■ eee 

Now we turn to /r —>■ eee decay. The branching ratio is given as [27, 59, 60] 


Br(/i ^ eee) = 


SGI 


1^ 


ND\ 


, , .n /'16, rrin 22\ 1. . 

+ \Ad? t'" - t + s ® 

' 3 rup 3 / 6 


i(2|F, 


L|2 


+ F, 


R\2\ 
Z \ ) 


-2A]y£)A}j + -AndB* 


+ h.c 


X Br(/i -)■ ez/gZ/^) (3.6) 

where Ad and And are the dipole and non-dipole contribution from the photonic 
penguin diagrams respectively. Also B represents the contribution from the box 
diagrams. Moreover, and are given as 

FzQl t?r FzQr 


pL _ ZiJL pi? _ 

^ Z 9 9*2/i Z 

sin 6w 


( 3 . 7 ) 


sin^0iv 

Here, Fz is the Z-penguin contribution and g^j^ and g^j^ are the Z-boson coupling to the 
LH and RH charged leptons respectively. In this model, Higgs penguin contribution 
will be suppressed by the small electron Yukawa coupling, and therefore we have 
only considered the photon penguin, Z-boson penguin and box diagrams. 

3.2.1 7 -penguin contribution 

First let us consider contributions from the photon penguin diagrams. In this case 
the 7 line of /uej vertex given in hgure 9 will have ee attached to it. The photonic 
non-dipole contribution. And for the doublet is in the following 


A doublet 
^ND 


3 


yeiVif^ f ^(n) / 


967r2 m 

i=l 


(3.8) 


The photonic non-dipole contribution, for the case of the quartet, will again have 
two parts, 

^quartet ^quartet . ^quartet /o q\ 

^ND — ^ND{n) ^ND{c) 

Here is the contribution of the neutral component and is the contri¬ 

bution of the charged component of the fermion triplet. 

( 3 . 10 ) 


/(quartet _ \ ^ \ ^ 2/eicr2/*/ro- 1 r^{n) (^ \ 

^ND(n)-2^ 2^ 9g^2 ^2^ 

(t=a7,AJ 


where again Xi^ = rnAo/rn?^. And the charged component of fermion triplet con- 

^ i 

tributes as follows, 

(3.11) 


^ quartet \ ^ \ ^ f pijc) / 

ND{c) 2-^ 2-^ 967r2 m? 

i=l a ^ 


with Xia = m‘t±/ml, and a = A++, S, A. The loop functions F^'^'>{x), F^'^\x), 

i 

GG) 

(x) and G^A(^x) are given in the appendix B. 
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3.2.2 Z-penguin contribution 

Now we focus on the Z-penguin diagram. The Feynman diagrams of one-loop contri¬ 
butions from the doublet and the quartet to the /reZ vertex are presented in hgure 
10. In Z-penguin diagram, the Z line of /reZ vertex will have ee line attached to 
it. For the doublet, the contribution is given by the neutral fermion. Following the 
formulas given in [60, 62, 63]^ 

1 ^ 

9ZC+C- C'24(mjv, , ruc) + me)] (3.12) 

i=l 


Here, gzc+c- is the Z boson coupling to (7^ of the doublet and g’‘j^ is the Z boson 
coupling to LH charged leptons given by 


9l = 


9 


1 

2 


-- -Fsin^6' 


w 


cos 9w 

On the other hand, the quartet contribution is 

^quartet _ ^quartet , r^quartet 

^Z{n) ^ ^Z{c) 

where the neutral fermion of the triplet contributes as 

^ Q 


(3.13) 


(3.14) 


T^quartet 

Z{n) 


Ibvr^ 


EE 22/eiail/i/4-72 9Zaia2 <724 (m^o, , m^j) 


2=1 <Ti,cr2 

+Vlin 


(3.15) 


where ai ^2 ^ and gzuia 2 is the Z boson coupling to cxi and cr 2 scalars of 

the quartet. The charged fermion of the triplet has the following contribution 

3 


T-iquartet 

^Z{c) 


Ibvr^ 


EE 

2=1 cri,o-2 


2/e2Cri 2/2 //(Ti 9 


2 C 24 {m^i,rnp±,Tnp±) -F - 


-F m^±(7o(m^j,mp,±,m^±) gzaia2 <724(m^±, m^^, 


+y*eia,y^^9^l9^Bl{mp±,m^^) 


(3.16) 


Here ai and a 2 range over the S, A, A++, and g^„±9^ is the coupling of Z boson to 

i i 

charged fermions. Moreover, Hi, Cq and C 24 are the loop functions, adopted from 
[60, 62, 63], presented in the appendix B. As Hi and C 24 arise from divergent loop 
integrals, for large M, 

1 1 

C 2 A{M,m,m) -)■ yln^ , Hi -)■ -In— (3.17) 

4 /i^ 2 /i^ 

^[60] contained a mistake in the calculation of Z-penguin diagram which was pointed out in [61]. 
Subsequently, correct results were presented in [62] and [63]. Moreover, Coo of [62] and C 24 of [63] 
only differ by an overall minus sign. 
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Therefore the combination 2 x02^ + yBi in Z-penguin contribution eq. (3.15) or in 
eq. (3.16) is vanishing at very large mass M when there are specihc relations set by 
group theoretical requirements in vertex factors x and y. 


3.2.3 Box contribution 

Lastly the box contribution for the doublet case, presented in hgure 11, is [60] 


^2 D doublet 


167r2 


E 

*5=1 . 


D. 


-y*eiyi^^y*ejyje + DomNimN^yl^yliyjf,yje 


(3.18) 


where, Dq = DQ{misf^,m]\f.,mc,mc) and Dq = DQ{mNi,mNj,fnc,mc) are loop func¬ 
tions given in the appendix B. 

For the quartet case, the contribution of the box diagram can be written as 


oquartet _ ^quartet , r^quartet 

“ ■“{**) ^ ■°(c) 


(3.19) 


with B 


quartet 

(n) 


is the contribution due to the neutral fermions and it is given by 


2 l/eiCTil/*A*<*'2l/ej(T2%'eo'l -^0ll^FPll^F9l/ei(Ti2/eio-2%’Mo'2l/jeo-i 

(3.20) 

where, Dq = DQ{mpo,mpo,m„^,m„^) and Dq = DQ{mpo,mpo,m„^,m„^). Here, cti _2 

’>‘3 ’^3 

ranges over and A^. 

The term corresponds to the contribution of the charged fermions and 

it reads 

BVr = E E (3.21) 

*5=1 0'1,0'2 

Here, Dq = DQ{mp±,mp±,ma^,ma^) and ui 2 ranges over A++, 5, A. 

3.3 fi — e conversion in nuclei 

The conversion rate, normalized by the muon capture rate is [27, 64-66] 


2 p quartet 

e ^(n) 


167r^ 


EE 

*5 = 1 0'1,0'2 


CR(/i — e. Nucleus) 


PeEemlGp 


^'em^eff^p 


Sti^ZT. 


capt 


\{Z + N){gpl + gpp) + {Z - N){g'iy + gpg)\ 


.(0) 


( 1 ) 


dl)M2 


\{Z + /V)(jS- + 4s) + (Z- Af)(4i( + 4>)|*} (3.22) 


Here, Z and N are the number of protons and neutrons in the nucleus, Z^jf is the 
effective atomic charge, Fp is the nuclear matrix element and Tcapt represents the 
total muon capture rate. Pe and Eg are the momentum and energy of the electron 
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(taken as ~ in the numerical evaluation), and Qx^j^ {X = 
in the above expression are given as 

9pK = 2 X] {9XKiq)Gx’^^ + gxK{q)Gx’'^^) 

q=u,d,s 

9xk = 2 ^ {.9XK{q)G%’^^ - gxK{q)G%''^^) 

q=u,d,s 


L, Ra.ndK = V, S) 


(3.23) 


gxK{q) are the couplings in the effective Lagrangian describing p — e conversion, 

Geff = {[gLSiq)eL^R + gRSiq)eRfJ,L]qq + [gLV{q)eLl ^+ gRV{q)eR'y^gR]q'lijiq] 

(3.24) 

are the numerical factors that arise when quark matrix elements are 
replaced by the nucleon matrix elements, 

{p\qTKq\p) = G^^'^^pTxP , {n\qTKq\n) = (3.25) 


For the inert scalar model, the g — e conversion rate receives the 7 , Z and Higgs 
penguin contributions. In 7 and Z penguin diagrams, qq (q=u,d,s) line is attached 
to 7 line of ge'j vertex and Z boson line of geZ vertex respectively. It doesn’t receive 
any box contribution because there is no coupling between inert scalars and quarks 
because of the Z 2 symmetry. Moreover, Higgs penguin contribution is small compared 
to 7 and Z penguin diagrams because of small Yukawa couplings thus neglected in 
our numerical analysis. The relevant effective coupling for the conversion in the inert 
scalar model is 


9LV(q) 9LV{q) 9LViq) 
gRV(q) = gLV(q)\L^R 
9LS(q) ~ 0 , gRS{q) ^ 0 


The relevant couplings are 


9LViq) ~ - Ad) 


9LV(q) 


V2gl + g1tFz 


G, 


nir 


(3.26) 

(3.27) 


Here Qq is the electric charge of the quarks and Z boson couplings to the quarks are 


9l - 


9 


-(T| - Qqsiv?ew) , Pr = 


9 


-Qq sin^ 9w 


f\ \ 6 y 1 jx f\ 

COS&'iy COsOw 

Also the relevant numerical factors for nucleon matrix elements are 

(Jfy — (jfy — Z , (_Ty — (_Ty — i 


(3.28) 


(3.29) 
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4 Results and Discussion 


In this section we have presented our numerical results and discussed the phenomeno¬ 
logical implications of those results for larger scalar multiplets. But before present¬ 
ing the results, we have listed all the constraints regarding dark matter and collider 
searches so that our analysis can focus on parameter space for where both inert 
doublet and quartet models are viable. 

There are two possible dark matter (DM) candidates in the inert scalar models. 
In the doublet model they are the lightest right handed neutrino, Ni and the lightest 
neutral scalar, S of the doublet. On the other hand, in the quartet model the neutral 
component of the lightest fermion triplet, and the lightest neutral scalar, S of the 
quartet can play the dark matter role. In both cases fermionic and scalar DM give 
rise to different phenomenology. In this preliminary study of comparing different 
LFV rates in inert scalar models, we have chosen the scalar as the DM particle and 
used the constraints associated with it in our analysis. 

4.1 Constraints and parameter space 

4.1.1 Collider constraints 

For the doublet scalar, the collider searches have put the following mass constraints, 
mc+ ~ 100 GeV, ms > 65 — 80 GeV and > 140 GeV [17-21]. Although there 
hasn’t been any collider studies on the quartet, one can recast the constraints of the 
doublet case onto the quartet. As the quartet scalar has the cascade decay channel, 
we can expect multilepton hnal states along with missing transverse energy similar 
to doublet. Therefore, the mass constraints for quartet, compatible with bounds on 
electroweak precision observable [67], are mAi 2 )^A++ ~ 100 GeV, ms > 65 — 80 GeV 
and mA > 140 GeV. Gonsidering S as the DM also set the mass hierarchy in quartet 
components: ms < rn^+ < m^++ < < mA- In contrast, the scalar masses in 

the TeV scale for both doublet and quartet are fairly unconstrained. 

In the case of fermions, the masses of RH neutrino in the doublet case are not 
constrained by current collider data. In contrast, fermion triplet of the quartet case, 
having gauge interaction, will have an accessible collider signature. In [68] the mass 
of the charged component of the triplet is excluded up to 270 GeV with 8 TeV 20.3 
fb“^ LHG data. Moreover, in [69] it was shown that the projected reach for 14 TeV 
collider with 3 ab“^ luminosity (High luminosity LHG phase) would he Mp < 500 
GeV, for (future) 100 TeV pp collider with 3 ab~^ luminosity in mono-jet searches, 
Mp ~ 1.3 TeV, and with 30 ab“^ luminosity, Mp < 1.7 TeV. 

4.1.2 DM Constraints 

The dark matter density of the universe measured by Planck collaboration is = 

0.1196 ± 0.0031 (68% GL) [70]. In the inert scalar model, there are two viable mass 
region of scalar DM. They are the low mass region {ms < my/) and the high mass 


13 


region (m^ 3> mw)- The low mass DM region of donblet model has been extensively 
stndied. In addition, same region for DM in the qnartet was addressed in [22] where 
it was shown that it is harder to achieve low mass dark matter with correct relic 
density compared to the donblet becanse, for most of the parameter space, bonnds 
on electroweak T parameter sets the mass of single charged component. A]*" close to 
the DM mass and therefore it is not only in tension with collider bounds but also 
opens up coannihilation channel and leads to a sub-dominant DM in the universe. 

In the high mass region of the doublet, as shown in [71], the DM mass starts from 
a lower bound of mo = 534 ± 25 GeV (where the thermal freeze-out only happens 
through the gauge interaction) to 20 TeV if the higgs-scalar coupling, \s < 27r. The 
maximal mass splitting compatible with correct relic density, are 

—'"i-sl ~ 16.9 GeV, |mc+— msl < 14.6 GeV (4.1) 

when ms 0(5 TeV). 

In the case of high mass region for the quartet, we have used FeynRules [72] to 
generate the model hies for MicrOMEGAS [73] and have found out that the DM mass 
starts from a lower bound of 2.46 TeV (freeze out only through gauge interaction)^ to 
upper bound of 14 TeV set again by \s ^ 27i bound. In this case, the mass splitting 
between the DM and other components are 

~ 16GeV, |m^+ — ms\ ^ 14GeV 
|mA++ — rns\ < 12GeV, |m^+ — ms\ ~ 1 GeV (4.2) 

when ms ~ 0(5TeV). hgure 2 presents the ms — plane with allowed region for 
both doublet and quartet scalar DM by the relic density and direct detection bound 
[74]. Here, Xs is effective coupling of S to Higgs held as can be seen in eq. (2.6). 
From hgure 2, we can see that there is an overlapping region on the plane where 
doublet and quartet DM satisfy the constraints simultaneously. 

The 7 coupling which controls the mass splitting between scalar (DM) and pseu¬ 
doscalar component, has the range 7 G [10“®, 2.7] for the doublet and 7 G [10“®, 1.36] 
to be consistent with the relic density. But it gets another constraint from bounds 
on DM inelastic scattering with nuclei. If the typical velocity of a DM particle, y is 
/3^c ~ 220 km/sec, the inelastic scattering is kinematically forbidden if the splitting 
between DM and the next to lightest component is larger, 

^ /S^m^iVfimcleus 

2 (m^ T Afnucieus) 

Therefore one would require, 7 > 10“^ to kinematically forbid the inelastic scattering 
of scalar DM with O(TeV) mass. As the inelastic scattering is mediated by the 

^without considering the Sommerfeld enhancement 
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exchange of Z boson and the scattering cross section is in the order of 10“^° — 
10 “^® cm^, which is mnch larger than the direct detection bonnds, the allowed range 
of 7 for doublet and quartet DM are 7 G [10“®, 2.7] and 7 G [10“®, 1.36], respectively. 




Figure 2. Correlation between the mass of the DM, ms and the effective coupling between 
the Higgs and the DM, A 5 for the doublet and quartet case. Here, the white region 
is excluded by the direct detection bound from the LUX collaboration [74]. The left 
figure represents the correlation without taking into account the Sommerfeld enhancement 
in the thermal freeze-out. In the right figure, for the green shaded region, Sommerfeld 
enhancement is not negligible. 


4.1.3 Gamma ray constraints and Sommerfeld enhancement 

Compared to the collider searches and DM direct detection experiments, indirect 
detection can set limits on the inert scalar DM at the TeV mass range because of a 
certain enhancement in the annihilation cross sections. 

At small relative velocity, two particles interacting via a long range force receive 
non-perturbative enhancement in the interaction cross section which is known as 
Sommerfeld enhancement [75]. When the mass of the DM is much larger than the 
mass of W and Z bosons, the electroweak interaction effectively behaves like a long 
range force, thus pair annihilation cross sections of the DM also receive Sommerfeld 
enhancements as pointed in [76-78]. At present, as the relative velocity of DM is 
about 10“^, Sommerfeld enhancement signihcantly boosts the indirect detection sig¬ 
nals, specially the gamma rays produced from the DM annihilation and put stringent 
constraint on the DM in the light of the experimental observations. In fact it was 
shown for the case of wino dark matter [81, 82] and minimal DM models (5-plet 
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fermion and 7-plet scalar with zero hypercharge) [83-85] (and references therein) 
that they are highly constrained to be the dominant DM of the universe by the ex¬ 
perimental limits on gamma ray spectrum due to the Sommerfeld enhancement in 
the pair annihilation cross section. 

Having electroweak charge, the heavy DM component of the inert scalar multiplet 
is also expected to have enhancements in both weak and scalar interactions. Although 
the full treatment of Sommerfeld enhancement for inert scalar model is beyond the 
scope of this work, following [79, 80], we introduce the dimensionless parameters to 
curve out the regions of the parameter space where the enhancement takes place 
and where the enhancement is negligible. The parameters are, = (udm/c)/q!, 
e<^ = (m<^/mDM)/tt and Here udm is the relative velocity of the 

DM particle, is the mass of the gauge boson carrying the force, 5 is the mass 
splitting between the DM and the next to lightest charged component of the multiplet 
and a is the coupling constant of the relevant interaction. It was shown in [80] that 
the Sommerfeld enhancement is relevant if e^, eg ^ 1. On the other hand, it is 
negligible for the region of parameter space where any of e^, > 1 . 

In the case of the minimal DM models, the processes contributing to the gamma 
spectrum from DM annihilation are, DM DM —>■ W^W~, ZZ where the decay and 
fragmentation of W and Z pairs produce secondary photons and DM DM —)■ 77 , 7 Z 
producing line spectrum of mono energetic photons. The Sommerfeld enhancement 
takes place when the DM-DM two particle state changes into DM^DM” two particle 
state, where DM^ is the next to lightest charged state, by exchanging W boson and 
subsequently charged states annihilate. For the minimal DM case, the DM and next 
to lightest charged state is almost degenerate (only loop induced mass splitting of 
the 0(100) MeV), so < 1 for = 1/30 and TeV scale DM and one can have 
Sommerfeld enhanced annihilation cross section. On the other hand, for the inert 
scalar models, the following terms in the scalar potential 

V D /3$V“<I)A'^T“A + 7 [($^er“$)(A^OT“A)^ + h.cj (4.3) 

can split the DM component and other charged component of the multiplet after 
electroweak symmetry breaking. For example, for quartet, when 7715 = 3 TeV and 
6 = m^+ — ms = 1.5 GeV, eg is 1.001. In addition, from figure 3, we can see that 
the bounds on electroweak precision observables allow maximum mass splitting to be 
8.78 GeV and corresponding eg is 2.46. Therefore for such mass splitting, according 
to [80], the Sommerfeld enhancement can be negligible in the inert scalar models. 

Moreover, Sommerfeld enhancement also affects the thermal freeze-out of the 
minimal DM as pointed out in [47, 48]. Such enhancement is also expected in the 
case of inert scalar DM. But if the thermal freeze-out happens after the electroweak 
phase transition, one can introduce enough mass splitting so that > 1. In fact, 
(5 = 1.5 GeV is compatible with the observed DM relic density of the universe with 
ms = 3 TeV for both doublet and quartet scalar DM. On the other hand, if freeze-out 
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S = m^ +—ms, GeV 


Figure 3. es vs 5 = m^+ — ms for DM mass, ms = 3000 GeV in the quartet. Here, 
blue points are allowed by stability conditions on the scalar potential and perturbative 
limits on scalar couplings. Red points are allowed by the bounds on electroweak precision 
observables. 

temperature, Tp is larger than the critical temperature of electroweak phase transi¬ 
tion, TpT, the thermal freeze-out takes place before the electroweak phase transition 
and there will not be any mass splitting to suppress the enhancement. Therefore 
thermal DM scenario of inert scalar DM will be different than that of the broken 
phase. But the value of the critical temperature of the electroweak phase transition 
depends on the model, order of the transition and its dynamics (see for example 
[86, 87]) . For this reason, we consider the range, Tpt = 100 — 200 GeV for the 
transition temperature. Now if Xf = Mdm/TV ~ 20, we see that for tudm > 4 TeV, 
freeze out takes place in the unbroken phase and will involve Sommerfeld enhanced 
annihilation cross-sections. 

On the other hand, for M^m < 4 TeV, the DM freezes out in the broken phase. 
So one can introduce the enough mass splitting, 5 ~ 1.5 GeV between the DM and 
next to lightest charged state to suppress the enhancement in the annihilation cross 
sections. 

For the inert scalar multiplets, apart from gauge interactions, the DM interacting 
via higgs exchange is also expected to have enhancement. In this case, the Yukawa 
potential I4c experienced by the DM is 

Vscir) = (Use - with (Use = (4.4) 

r 4:71 Trig 

For example, if ms = 1 TeV and Xs = vr, Osc = 0.047 so = 2.6 for the Higgs 
exchange, therefore the enhancement is generally not important for scalar interaction 
with the DM mass at TeV range. 

In summary, although the DM with mass at TeV range in the inert scalar model 
is expected to have Sommerfeld enhancement in the gauge interactions and can have 
signihcantly enhanced indirect detection signal, there is a small common region of 
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parameter space for doublet and quartet as seen from figure 2 (right) where one can 
have enough mass splitting to suppress the Sommerfeld enhancement in the inert 
scalar models and such mass splitting is compatible with the observed DM relic 
density. Therefore in the subsequent analysis, we only focus that small region of 
parameter space with benchmark point, ms = 3 TeV and S = 1.5 GeV and have 
left the complete analysis of Sommerfeld enhancement in the inert quartet case for 
future work [90]. 

4.1.4 Scalar coupling and LFV rates with scalar DM 

There is a correlation between the 7 coupling of the scalar sector and the rate of 
LFV processes when R in eq. (2.1) is a real orthogonal matrix. As we can see from 
eq. (2.15) and eq. (2.17) that the smaller value of 7 leads to smaller value of the loop 
factor Aj and thus neutrino mass. This in turn increases the Yukawa coupling, as in 
eq. ( 2 . 20 ), and becomes inconsistent with perturbativity bound eq. ( 2 . 21 ) when 7 is 
very small. On the other hand, large value of 7 implies larger separation in m 5 and 

and also in m^+ and m^+, thus larger value of Aj and in this case the value of 
Yukawa coupling is reduced. In hgure 4 (left). We have illustrated this by comparing 
Br(/i —>■ 67 ) for 7 = 10“® and 10“® respectively. We can see that for 7 = 10“^, the 
rate has become out of reach for current and future experiments. Therefore in the 
case of real R matrix, 7 ~ 0(10“®) leads to appreciable LFV rates. However we 
have seen in Sec. 4.1.2 that as one would require, 7 > 10“® to kinematically forbid 
the inelastic scattering of scalar DM with O(TeV) mass so considering only real R 
will lead to negligible rates of LFV processes. 

On the other hand, in the case of complex R, such correlation between 7 and the 
rates of LFV processes is not straightforward because the size of Yukawa coupling 
also depends on the imaginary part of the complex angles in R. For simplicity, we 
have added an imaginary part, Im(z), in three angles of R and in hgure 4 (right), 
we can see that, despite having 7 = 10“^, Br(/r —)■ 07 ) become comparable to the 
current bound with increasing values of Im(z). Again perturbativity of the Yukawa 
coupling typically put upper bound on Im(z) of 0(3 — 5). Therefore, one can have 
viable scalar DM in both doublet and quartet models where ^ > 1 with appreciable 
LFV rates by tuning the Imaginary part of complex angles in R. 

4.1.5 Viable parameter space 

The parameter space for the model consists of {Mq, a,/3, 7 } of the scalar sector and 
of the fermionic sector. Here Mat and Mp are the masses of RH neutrino 
and real fermion triplet (as the components of the triplet are degenerate at tree level) 
respectively. 

The focus of this preliminary study is the comparison among different LFV rates 
in doublet and quartet model with scalar DM. At hrst, from hgure 2, as an exemplary 
point, we have chosen the mass of scalar DM to be ms = 3 TeV with A 5 = 1.3 in 
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Normal Hierarchy 


Normal Hierarchy 



Figure 4. Left figure presents the dependence of the rate of LFV processes on the 7 when 
the Risa real matrix. Here we have considered only Br(/i —?■ ey) for illustration. The brown 
and blue represents the rate in the doublet and quartet cases respectively for 7 = 10“®. On 
the other hand, the orange and red points represents the rate in the doublet and quartet 
cases respectively for 7 = 10“^. Right hgure presents the correlation of the rate in doublet 
(brown points) and quartet (blue points) with imaginary part of the complex angle, Im(z), 
when we consider complex R matrix. Here The scalar mass is fixed at m-scaiar = 3000 GeV 
and 7 = 10“®. The black horizontal line is the current bound 5.7 x 10“^^ and red line is 
projected bound 6 x 10 “^^. 

the ms — Ag plane so that scalar DM is viable both in donblet and qnartet model. 
Moreover, 7 is set to be 10“® to be consistent with bounds from DM direct detection. 
As the components of the scalar multiplet are almost degenerate apart from the very 
small splitting induced by non zero 7 . Therefore we set the average mass of the 
scalar components at ruscaiar = 3 TeV. 

There are two sets of fermion mass range we have considered in our analy¬ 
sis. For the comparison of LFV rates with the variation of fermion masses both in 
doublet and quartet model, we have evaluate them in two sets, namely, i) where 
^ < 1 so that the scalar component ceases to be the DM and ii) 

where > 1 where the scalar component is the DM. We have varied the masses of 
RH neutrinos and the fermion triplet within the range, Mjv(f) ^ (270 GeV, 30 TeV) 
which encompasses both sets mentioned above. 270 GeV is taken as the lower limit 
of fermion mass as triplet fermion is excluded up to that mass in collider searches. 
Also such range is considered to see how the LFV rates vary with the mass of the 
fermion in addition to the DM aspects of inert scalar model. 

We have used the experimental values of low energy neutrino parameters, VpmnS) 
as the input in eq. (2.1) for Yukawa couplings. For both normal 
and inverted hierarchies, we could only vary the lowest neutrino mass, , the Dirac 
phase, 5 and Majorana phases, and three complex angles, Zi, Z 21 2^3 of, R. In 

our numerical analysis,as an simplihcation, the lowest neutrino mass is set to = 1 
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meV, 6 G [ 0 , 2n], a = (3 = 0 and common imaginary part in Zi = 6i + i Im(z) 

with the range (0, 5). 

Snmmarizing, our input parameters in the numerical scans are {Mq, a, (3, 7 , Mjsf = 
Mp = M, rui/j, 6, 6 * 2 , 6 * 3 , Im(z)} satisfying all the constraints mentioned above 

Therefore, we can compare the LFV rates in both models for common viable point 
in the parameter space. 

4.2 LFV processes 

In the inert scalar models with scalar DM in the high mass regime, there is no direct 
correlation between the Yukawa couplings and DM properties. Also we have seen 
that the real matrix R and 7 ^ 10“^ (scalar DM direct detection constraint) give 
rise to small Yukawa couplings which in turn lead to LFV rates beyond the reach 
of current and future experiments as seen in hgure 4 (left). But the size of the 
Yukawa coupling can be enhanced by varying the imaginary part of complex angles 
of R without substantially affecting the phenomenology of the scalar DM and despite 
having 7 > 10“®, we can easily obtain the LFV rates within the experimental range. 

So hrst we have compared the rates of /r —>■ ey, /i —>• eee and /i —e conversion rate 
with 7 = 10“® and the real R matrix by varying the fermion masses for the doublet 
and quartet models. Then we vary Im(z) within its constrained limits and determine 
the region allowed by current and future bounds on the rates of these three LFV 
processes for both doublet and quartet cases. Also when ^ > 1 , we have our scalar 
dark matter in both doublet and quartet models. 

4.2.1 Br(/r —)■ ey) 

Due to the excellent bound put by the MEG collaboration [33, 34], /r —>■ ey is one of 
the most well studied LFV processes, hgure 5 shows the comparison of this process 
between the doublet (brown points) and the quartet (blue points) scalar. We can 
see that the quartet contribution to /i —)■ ey is larger than that of the doublet. For 
the same parameter point, in the quartet case, additional charged and neutral scalar 
(A]*", Af, A^^, S and A) and fermion states (F° and F^) enter in the loop compared 
to single charged scalar (G^) and neutral fermion state (Yj) in the doublet case and 
as the contributions of extra states are additive, the rate has increased in the quartet 
case than that of the doublet. From hgure 5 we can see that Br(z/ —)■ ey) is larger 
for the quartet than the doublet for both < 1 and ^ > 1 (where the doublet and 
quartet scalars are the DM). 

4.2.2 Br(/i —)■ eee) 

In /r —)■ eee, the dominant contributions are coming from y-penguin and Box di¬ 
agrams. The Higgs penguin diagram is suppressed by the small electron Yukawa 
coupling. The Z penguin contribution is small because of the cancellation that takes 
place between G 24 and Bi terms in eq. (3.12) and also between the same terms in eq. 
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Normal Hierarchy 


Inverted Hierarchy 



Figure 5. Correlation between ^ = -^^(p)/?7iscaiar Br(^ — )■ 67 ) for doublet (brown 
points) and quartet (blue points) with normal (left fig.) and inverted (right fig.) hierarchy 
for light neutrino mass. Here we have taken to be degenerate, random Dirac phase 

6 and random real matrix R. Also we have set Majorana phases Oiy and /Jj/ to be zero in 
this case. The scalar mass is fixed at mscaiar = 3000 GeV. Also 7 = 10“® and light neutrino 
mass, m-i/j = 1 meV. 

(3.15) when ma-^ = ma- 2 - Moreover, similar cancellation takes place between the hrst 
two lines and third line of eq. (3.16) due to the specihc relations among the couplings 
in front of the vertices. Therefore Z penguin contribution is also small in /i —)■ eee 
for both inert doublet and quartet case. Also note that the Z contribution in the 
quartet case is relatively bigger than that in the doublet because in the quartet 
and ma -2 exactly equal when ui 7 ^ a 2 . Hence one receives larger Z-penguin 

contribution in the quartet compared to the doublet. Still this contribution is nu¬ 
merically not signihcant if we compare it with 7 penguin diagram or box diagram 
contributions. From hgure 6 , we can see that Br(p —>■ eee) is larger for quartet (blue 
points) compared to the doublet (brown points) for both ^ < 1 and ^ > 1 cases. 

4.2.3 fi — e conversion rate 

Another prominent LFV process currently under investigation is the /i — e conversion 
in nuclei. Here we have calculated the /i — e conversion rate for Ti and An nuclei 
in the inert model with doublet and quartet. From hgure 7, we can see that the 
/i — e conversion rate is larger for the quartet (blue points) compared to the doublet 
(brown points). The dip occurs in the doublet contribution at ^ = 1 because at that 
value, the dipole contribution and the non-dipole contribution are 

equal as they are coming from single 7 penguin diagram involving charged scalar 
and neutral fermion Ni and eq. (3.26) indicates that the effective coupling is zero 
for doublet at that point. On the other hand, for quartet case and at 

.^ = 1 are different because more than one charged scalar contribute to the 7 penguin 
diagrams. Again we can see from hgure 7 that the conversion rate is larger for the 
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Figure 6. Correlation between | /iTT-gcaiar Br(// —)■ eee) for doublet (brown 

points) and quartet (blue points) with normal (left fig.) and inverted (right fig.) hierarchy 
for light neutrino mass. Here we have taken same input parameters as in Br(/i —?• ey). 


quartet than that for the doublet for both < 1 and ^ > 1 cases. We have not 
included the hgure for fi — e conversion rate in An nuclei as it is similar to hgure 7. 


Normal Hierarchy 


Inverted Hierarchy 



^ — M. N I Scalar 


^ — M. N I ffi Scalar 


Figure 7. Correlation between ^ and ^ — e conversion rate for Ti nucleus, 

for doublet (brown points) and quartet (blue points) with normal (left fig.) and inverted 
(right fig.) hierarchy for light neutrino mass. Here we have taken same input parameters 
as in Br(/r —)• ey). 


4.2.4 LFV rates in the doublet and quartet 

As expected, the LFV rates seen in hgure 5, 6 and 7 are very small for real R and 
7 = 10“^. The rates will reduce even more if we increase 7 . Still the rates are larger 
for the quartet compared to the doublet for .^ < 1 and ^ > 1 case where scalar is 
treated as the DM candidate. Now we increase the value of Im(z) and calculate the 
LFV rates with increasing values of M. 

From hgure 8 , we can see that LFV rates in the quartet are more constrained 
than those in the doublet for common parameter space satisfying all the restrictions 
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Figure 8. The Im(z) plane for degenerate random Dirac phase 6 , zero Majorana 

phases Uu = = 0 and light neutrino mass, = 1 meV. The scalar mass is m-scaiar = 

3000 GeV with 7 = 10“®. In (left), the current bounds are imposed: Br(/r —)■ e'y) ^ 
5.9x 10“^^, Br(/x —>■ eee) ^ 1 x 10“^^ and /r—e conversion rate for Ti < 4.3x 10“^^. In (right), 
the future sensitivity are considered: Br(/r —)• ey) 6.4 x 10“^^, Br(/i —>■ eee) 1 x 10“^® 
and n — e conversion rate for Ti < 10“^®. 

of Sec. 4.1. The allowed regions on ^-Im(z) plane are reduced further for both 
doublet and quartet models if one imposes the sensitivity of future lepton flavor 
violating experiments. The case for inverted hierarchy shows similar pattern so we 
have only presented results regarding normal hierarchy. 

5 Conclusions 

The scotogenic model is a well studied neutrino mass model and lepton flavor vio¬ 
lation is one of its important phenomenological aspects. In this study we present 
the comparison among different LFV processes in the inert doublet and the quartet 
model, taking into account the current experimental limits and future sensitivity. 
There are two possible dark matter candidates in the inert scalar models: scalar 
and fermionic DM. In this study we have considered scalar DM and evaluated LFV 
rates for common parameter space subjected to collider bounds, DM constraints for 
doublet and quartet model and low-energy neutrino parameters. Our results are 
summarized as follows 

• Br(/i —)■ ey), Br(/r —)■ eee) and fi — e conversion rates in nuclei in the quar¬ 
tet model are larger than those in the doublet model for the same parameter 
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space as seen from figure 5, 6 and 7. In the case of higher scalar representa¬ 
tion more particles enter into the loops and their contributions are additive in 
the LFV processes. Therefore we can have larger rates of different LFV pro¬ 
cesses compared to the lower scalar representation. From hgure 8, we can see 
that, LFV processes in higher scalar representation are more constrained by 
the current and near-future experiments. In addition, this phenomenological 
result is complementary to the appearance of low scale Landau pole for higher 
representations [22, 88, 89]. 

• There is no signihcant deviation from hgure 5-8 for non-degenerate right handed 
neutrinos and real fermion triplets. In the case of large hierarchy, mTVg >> 
mNi 2 J dominant contribution comes from only the lightest generation. 

We would like to emphasize here that the conclusion of our preliminary study 
is applicable to the inert scalar models where scalar DM is considered. But there is 
much room for an improved analysis. For example, in the case of fermionic DM, the 
DM constraints will be different and will have different viable parameter set for the 
LFV rate comparison. Also one needs to study the DM properties and viability of a 
common parameter space where ^ ~ 1. Therefore, further quantitative analysis of the 
fermionic DM aspects in the quartet model will be presented in a future publication 
[90]. Furthermore, similar analysis can be carried out for r —?• /xy, r —>■ eee, r —?• /i/i/I 
in the inert scalar models to probe the flavor structure of the Yukawa sector and to 
have better constraints on the higher scalar representation in the light of experimental 
limits. 
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A Scalar masses 


A.l Inert Doublet 

The mass spectrum for the inert doublet in our parametrization eq. (2.2) is, 

^ + i/3 + 7^ 

m\ = + i + i/3 - 7 ^ 

m^Q = Mg + 2 (^-1) 


A.2 Inert Quartet 

In the inert quartet case, the 7 term, apart from splitting S and A, also mixes two 
single charged components of the quartet. According to eq. (2.9), the mass matrix 
for single charged helds in (A"*", A'"*") basis is 





Mq +l{a + \(5)v‘^ 


(A. 2 ) 


Diagonalizing the mass matrix, we have mass eigenstates for single charged helds, 
A/^ = A+ cos 9 + A'+ sin 6*, A^ = — A+ sin 9 + A'+ cos 9 with tan 29 = — _ 

Therefore the mass spectrum of the quartet is 


^S(A) = + 27^ 

m\++ = M^ + ^(^a- i/3^ 

^l+(A+) = ^0 + ^ T i\//^2 + i272^ (A.3) 


Because of the mixing between two single charged states, the mass relation is 


2,2 2,2 

rus + = m^+ + m^+ 


(A.4) 
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B Loop functions 


The loop functions relevant for the dipole and non-dipole form factors from f^e'y 
vertex are 


F^^\x) 

G^^\x) 

G^^\x) 


1 — 6a; -I- 3a;^ -|- 2x^ — 6a;^lnx 

6(1 -a;)4 

2 -|- 3a; — 6a;^ -|- a;^ -|- 6a;lna; 

6(1 -a;)4 

2 — 9a; -|- 18a;^ — lla;^ -|- 6a;^lnx 
6(1 -a;)4 

16 — 45a; -|- 36a;^ — 7x^ + 6(2 — 3a;)lna; 
6(1 -a;)4 


In the following we collect the Passarino-Veltman loop functions. 




1 m\-ml + 2m\ln^ i ^2 

2 4(mf — 2^"^ fF 


(B.l) 

(B.2) 

(B.3) 

(B.4) 


(B.5) 


Go{mi,m2,m3) 


2 / 2 2 M 7ni 

m2\mi — 7723 ) 111 ^ — 


{m{ — m\){m\ 


{m\ — 

ml){ml - ml) 


(B.6) 


G2i{rni,m2,mFj 




—2{m\ -I- m\)ml 



1 


m 


1 


{ml - ml) 
(B.7) 


Do{mi,m2,m3,m4) 



{ml — ml){ml — ml) {ml 



{ml — ml){ml — ml) {ml 





{ml — ml){ml — ml) {ml — ml) 


(B.8) 


Ilo(mi, 7712, m 3 , 7714 ) 



{ml — ml){ml — ml) {ml 



{ml — ml){ml — ml) {ml 





{ml — ml) {ml — ml) {ml — ml) 


(B.9) 
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C /ie7 vertex, ^leZ vertex and box diagrams 

C.l /ie 7 vertex 

Here we present in figure 9 the Feynman diagrams of one-loop contributions of the 
doublet and quartet to the /rey vertex. 



Figure 9. pey vertex and the self energy diagrams of the external fermions for the doublet 
(first row) and the quartet cases (second to fourth rows). 


C.2 jjieZ vertex 

We present in figure 10 the Feynman diagrams of one-loop contributions of the dou¬ 
blet and the quartet to the /reZ vertex. 
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Figure 10. [leZ vertex and the self energy diagrams of the external fermions for the 
doublet (first row) and the quartet cases (second to fourth rows). 

C.3 Box diagrams 

The box diagrams for the doublet and the quartet cases are given in hgure 11, 



Figure 11. Box diagrams for the doublet (hrst row) and the quartet (second and third 
rows). 
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